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ABSTRACT 

Context. 

Aims. The main goal of this paper is to prove that accreting very low-mass stars and brown dwarfs can be identified in IPHAS, a 
Hcf emission survey of the northern Milky Way.Full exploitation of the IPHAS database and a future extension of it in the southern 
hemisphere will be useful to identify very low-mass accreting objects near and far well-known star forming regions. 
Methods. We have used Virtual Observatory tools to cross-match the IPHAS catalogue with the 2MASS catalogue. We defined 
photometric criteria to identify Ha emission sources with near-infrared colours similar to known young very low-mass stars and 
brown dwarfs. 4000 candidates were identified that met our criteria over an area of 1600 square degrees. We present low-resolution 
optical spectra of 113 candidates. Spectral types have been derived for the 33 candidates that have spectroscopically confirmed Ho- 
emission, negligible reddening and M spectral class. We have also measured Ha emission and investigated the Nal doublet (818.3 
nm, 819.5 nm) in these 33 objects. 

Results. We confirm that 33 IPHAS candidates have strong Ha indicative of disk accretion for their spectral type. 23 of them have 
spectral class M4 or later, of which 10 have classes in the range M5.5~M7.0 and thus could be very young brown dwarfs. Also many 
objects have weak Nal doublet, an indication of low surface gravity. 

Conclusions. We conclude that IPHAS provides a very valuable database to identify accreting very low-mass stars and brown dwarfs, 
and that Virtual Observatory tools provide an efficient method for identifying these objects over large areas of the sky. Based on our 
success rate of 23 Ha emission objects with spectral type in the range M4-M7 out of 113 candidates with spectroscopic follow-up, 
we estimate that there could be hundreds of such objects in the full IPHAS survey. 

Key words, very low mass stars 



1. Introduction 

Since the first unambiguous discovery of b rown dwarfs (BDs) 
(iRebolo et alJ Il995t iNakajima etlD 1 19951) this field has pro- 
gressed rapidly. Considerable observational effort has been de- 
voted to identifying BDs in the known nearby star-forming re- 
gions (SFRs) and young open clusters. These objects provide 
crucial information on the dependence of key physical proper- 
ties, such as disk properties, multiplicity and the shape of the 
Initial Mass Function on primary mass. 
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Due to mass accretion processes, many young low-mass 
stars and BDs show Ha emission stronger than the emission 
expected from chromospheric activity. Studying the Ha equiv- 
alent width and the spectral type using low-resolution spec- 
tra, it can be determined whether o r not a star is accreting 
(iBarrado y Navascues & Martini 120031) . Thus Ha surveys have 
the potential to identify very young stars and BDs that are still 
accreting from their disks. 

Some Ha searches for very low-mass (VLM) stars and BDs 
have already been instrumental in detecting these objects, such 
as for example those carried out in high latitude molecular 
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cloud s (lMartin&KuDlll996h or in Orion OBI dBriceno et alj 
1200 lb . Evidence of accretion disks and infrared excess in 
young VLM stars and BDs has been confirmed in several 
SFRs, such as p Op hiuchus ( Wilking et al. 1999), the Trapezium 
Muench et al.l l2001.) . Chamaleon I CNatta & Testi 2001), IC348 



( Java wardha na et al] 20031) or cr Orionis (Oliveiraet al.i 120021: 
IZapatero Osorio et alJ l2002 a': Muz erolle et al.i |2003). Some 
spectroscopic studies have emphasised on the analysis of the Ha 
emission of substellar objects a nd led to the detection of v ery 
cold obje cts such as S Ori 55 (|Zapatero Osorio et al.ll2002bl) or 
S Ori 71 jBarrado v Navascues et al.ll2002l) . 

The INT Photometric Ha survey of the Northern Galactic 
Plane (IPHAS) is a valuable source for discovering young VLM 
stars and BDs using Ha emission. It covers 1800 square degrees 
of the northern Milky Way in the latitude range -5°< b < 5°. It 
provides (Sloan) r', i' and narrowband Ha photometry down to 
a magnitude limit of r'~ = 20 (10 cr). The data is taken using the 
Wide Field Camera (WFC) on the 2.5m Isaac Newton T elescope 
(INT) (iDrew et al.ll2005l:lGonzalez-Solares et al.ll2007h . 

So far the overwhelming majority of the surveys for young 
VLM objects are concentrated in the known SFRs and nearby 
young clusters.Tak e as an example recent resea rch in Taurus 
dGuieu et all |2006|) . Or ion OBI A and OBIB jPownes et alj 
I200d) . Cha maeleon I jLuhmanl l2004al) . Cha maeleon II and 
Ophiu chus dAllers et al.i l2007h or crOrionis dCaballero et alj 
I2007h . 

The IPHAS survey offers a complementary approach 
because it allows us to use Ha as a primary selection criterion, 
and it provides a wide area coverage around or outside the 
well-known SFRs and clusters. 

The Virtual Observatory (VO) is a recent initiative with the 
goal of managing large databases in an organised manner in or- 
der to make an efficient use of astronomical archives. In this 
work VO tools have been applied for the first time to the search 
for young VLM objects. We have made a cross-correlation of 
an early version of the IPHAS point source catalogue with the 
Two Micron All Sky Survey (2MASS) point source catalogue 
jSkrutskie et al. 2006). The rest of the paper is organised as fol- 
lows: Section 2 describes the selection criteria used for identifi- 
cation of young VLM candidates. Section 3 deals with spectro- 
scopic follow-up observations of 113 candidates selected using 
our criteria. Section 4 discusses the analysis of the spectra lead- 
ing to the determination of spectral types and equivalent widths. 
Section 5 presents the main results of our study. 



2. Sample selection 

The VO offers the possibility of efficiently cross correlating 
large multi-wavelength databases. AZaii/nQis an interactive soft- 
ware sky atlas allowing the user to visualise digitised astronomi- 
cal images, superimpose entries from astronomical catalogues or 
databases, and interactively access related data and information 
from the Simbad database, the VizieR service and other archives 
for all known sources in the field. 

We have used AZflii/« to look for new VLM objects via cross- 
correlation of the IPHAS catalogue with 2MASS. The photomet- 
ric data available in these catalogues provides magnitudes in Ha, 
r', i', J, H and K for the selected objects. 

We used an early version of the IPHAS point source catalogue 
that covered about 1200 square degrees. Our search was re- 
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stricted to the RA range from 18 to 05 hours. We identified 4000 
candidates that met the following criteria: 

- IPHAS-2MASS coincidence in coordinates within 1" . 

- The IPHAS sources should be classified as stellar or prob- 
ably stellar in the r', i' and Ha bands and have colours 
in the range 1.1 < r'-Ha < 3.0 (to avoid potential arti- 
facts and s e lect o bjects with Ha in emission according to 
iDrew et all (|2005|) who show that M6 dwarfs have IPHAS 
colours r'-Ha = 1.06 for E(B-V) = 0.0, and r-Ha =1.14 
forE(B-V) = 1.0); and magnitudes i' < 18.5. 

- The 2MASS sources should have qflag A or B; colours in the 
range 0.7 < J-H < 1.3 (to discard strong reddened objects 
and red giants); and 0.4 < H-K < 1.1 (for selecting cool 
photospheres and/or infrared excess). 



3. Follow-up spectroscopic observations 

We have carried out 3 campaigns of follow-up low-resolution 
spectroscopic observation of our targets. Altogether, spectra for 
113 candidates have been obtained, which represents only 3% 
of our total sample. We now discuss these runs in chronological 
order: 



3.1. William Herschel Telescope observations 

Our first run took place on August 1-2, 2006 with the ISIS long 
slit spectrograph on the 4.2m William Herschel Telescope in La 
Palma as part of the international time program led by Janet 
Drew. The R158R grating in red arm was used. The instrumen- 
tal setup gave a dispersion of 1.63 A/pixel with a wavelength 
range of 5400-10300A. We took spectra with exposure times 
ranging from 700 to 1800s for 35 young VLM candidates. Due 
to poor weather conditions in the first night (variable seeing of 2- 
3") the slit width was 2.0"(FWHM ~ 15A). In the second night 
we had much better seeing, ranging from 0.7"to 1", and hence 
a slit width of 1.0" was used (FWHM ~ 7A). Each night two 
flux calibration standard stars were observed, and Ne-Ar arcs 
and lamp flat-fields were taken at the beginning of the nights. 
Bias and flat-field correction has been applied to all our CCD 
frames. Wavelength calibration and instrument response correc- 
tions were made on the science spectra using standard tasks in 
the IRAF environment. 



3.2. Northern Optical Telescope observations 

Our second run occurred on October 10-14, 2006 with the 
ALFOSC spectrograph on the 2.5m Nordic Optical Telescope in 
La Palma. Grism 5 was used corresponding to a dispersion of 3. 1 
A/pixel and a wavelength range of 5000-10250A. A slit width 
of 1.0"(FWHM ~ 16A) was used and exposure times ranged 
from 300 s to 2400 s. 56 new candidates were observed. One 
flat-field and one arc lamp exposure were taken immediately be- 
fore or after each science target. A flux calibration standard star 
was observed at the beginning of the nights. Bias correction was 
performed using an average bias for each night, and flat-field 
correction was made with the normalised flat-field exposure ob- 
tained at the same telescope position as the scientific spectrum. 
This procedure improved the correction of the fringing pattern in 
the CCD detector response over the correction that was obtained 
using an average flat-field. The calibration in wavelength was 
done using the images of He-Ne-Th-Ar arcs obtained after each 
scientific image. Tests were done calibrating with the emission 
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lines of sky that appear in the spectra instead of with the arcs, 
nevertheless the results of the calibration were not better. Tests 
were also made using different polynomial functions to fit the 
wavelength solution. Finally, a 4th order Legendre polynomial 
was chosen. We measured that the displacements of known sky 
airglow emission lines with respect to the expected positions are 
generally less than ~2A , which is acceptable because it is less 
than 1/8 of our spectral resolution. The standard HD227900 was 
used for flux calibrations. All data reductions were carried out 
using IRAF. 

3.3. Lick Observatory Shane Telescope observations 

In July 2007 we used the KAST spectrograph on the 3m Shane 
Telescope at Lick observatory with the 600/7500 (2.32 A/pixel) 
and 300/7500(4.6 A/pixel) gratings in the red arm only with a 2" 
slit width and exposure times from 300s to 2400s for a total of 25 
candidates. The wavelength range covered was of 5600-8500 A 
for the 600/7500 grating and 5000-10500 A for the 300/7500 
grating. Ne arcs and flat-field lamp exposures were taken for 
each position of the scientific images. Part of the observations 
of the last night were affected by adverse weather (cirrus). In 
general the quality of the data collected during this campaign is 
poorer than the two previous runs. During the first two nights 
the 600/7500 grating was used (FWHM ~ 5A). Although the 
weather was good, with seeing of ~1.5", the spectra were very 
noisy and weak, sometimes so much that it was not possible to 
extract them. For the 2 last nights we decided to change to the 
300/7500 grating (FWHM ~ 9A) in order to improve the sig- 
nal to noise ratio of the spectra. For the reduction, flat-field nor- 
malised images obtained after each scientific exposure were ap- 
plied. For wavelength calibration Ne arcs also taken after each 
exposure was used. Depending on the night and the grating, dif- 
ferent spectrophotometric standards were used for flux calibra- 
tions. All data reductions were carried out using IRAF. 



4. Spectral analysis 

We found that 42 of our 113 targets have strong Ha emission. 
Therefore we have had a 37% success rate in the confirmation 
of Ha emission among our candidates. We find 12 of the 42 con- 
fir med targets in the IP HAS catalogue of emission-line sources 
of lWitham et alj (|2008). For the sample selection we used a pre- 
liminary version of the IPHAS catalogue which was available to 
us in July 2006. Due to the improvement of the Ha zero point 
calibration on the IPHAS catalogue, although all of our candi- 
dates where selected as likely Ha emitters with r'-Ha > 1.1, 
most of the non Ha emitte rs are out of the ma in stellar locus 
of emission-line sources of lWitham etall ( l2008i) in the IPHAS 
colour-colour plane, and only 48 of them fit the selection criteria 
according to the new photometric IPHAS catalogue. According 
to this, we have 48 observed candidates that really fit the crite- 
ria and from those 40 of them show Ha emission which is 83% 
success rate. In the following analysis we only retain the objects 
with confirmed Ha emission. 

4.1. Spectral type determination 

The spectral type has been determined using 3 different methods: 
(a) the Hammer code; (b) fitting with template spectra; and 
(c) measurement of spectral indices that have a known relation 
to spectral classification. 
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Fig. 1. An example of spectral type determination using the 
method. The best fit (M6) is represented in red colour overlaid on 
the observed spectrum J210404.87-I-535124.4 (black). For com- 
parison we represent also the M5 (blue) and M7 (green) refer- 
ence objects. The spectra are normalised at 8123A. 



4.1 .1 . Comparison with template spectra: Tine Hammer code 

Ha mmer^ is a tool wri tten in IDL for the analysis of stellar spec- 
ti-a (Covey et al."2007'). With this routine each spectrum is anal- 
ysed automatically to predict the spectral type by calculating a 
list of spectral indices (measuring the strength of CaH and TiO 
bands) to a library, considering the uncertainty of each index 
and also do measurements of Ha as tracer of activity in late type 
stars. In addition, it allows the user to manually compare and to 
change the final spectral type that is initially assigned automati- 
cally by the program via a rough comparison with the sample of 
reference spectra. The routine has been applied to all our targets 
with confirmed emission in Ha. In most cases, it was necessary 
to select interactively the spectral type. Nevertheless, applying 
this method is useful because it allows to quickly inspect differ- 
ent comparisons of the targets with templates and to derive a first 
estimate of the spectral type. 

A.\.2.x^ metliod 

An IDL program has been written to calculate the value com- 
paring our target spectra with those of M-type reference field 
dwarfs. We found the best spectral type by minimisation of the 
X^ value. As reference for the spectral classification, stars of th e 
I AC ultracool dwarf catalogue were used (Martin e t alj|2005h . 
We selected one dwarf for each spectral subclass between MO 
and M9. The dwarfs selected were: SDSS J125350.99-001010.3 
(MO), SDSS J004830.1 1-003029.1 (Ml), SDSS J000201.55- 
010636.5 (M2), SDSS JOl 11 13.97-003128.4 (M3), SDSS 
J002328. 15+001709.02 (M4), SDSS Jl 13922.08+000048.8 
(M5), SDSS J150026.34-003928.0 (M6), SDSS J100218.61- 
000408.46 (M7), SDSS J135859.02-005357.5 (M8), and SDSS 
J024958. 36-003409.99 (M9). The spectra of the catalogue 
were resampled using the task linterp and normalised in the 
wavelength range of /l=81 10-8136 A (central wavelength 
/1=8123 A) to make the calculation of x^- 



^ See |http://www.cfa.harvard.edu~kcovey/thehammer| 
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Table 2. Spectroscopic results. 
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IPHAS name/position 


SpT 


SpT 


SpT 


SpT 


W(Ha) 




J[RA(2000)+Dec.(2000)] 


{Hammer) 


Or") 


(indices) 


(adopted) 


(Hammer) 


W(Hq') 


J001649.57+654241.8 


M5 


M7 


M4.5 


M5.5±1.0 


140.3 


250±80 


J001655.91+654732.8 


M6 


M7 


M5.5 


M6+0.5 


57.7 


74±15 




M5 


M7 


M5.5 


M6+1.0 


96.1 


133±15 


JOl 1443.00+620820.9 


M4 


M5 


M4.5 


M4.5+0.5 


192.7 


276±50 


J012348.67+614931.8 


M4 


M4 


M4.5 


M4±0.5 


135.6 


230±30 


J013720.01+645957.7 


M4 


M4 


M4.5 


M4±0.5 


112.7 


170±30 


J023616.00+615609.6 


Ml 


M4 


- 


M2.5+1.0 


133.7 


205+30 


J035449. 17+530903.3 


M3 


M4 


M2.5 


M3+0.5 


93.7 


164+30 


J035823.95+522312.6 


Ml 


M2 


- 


Ml.5+0.5 


147.2 


214+30 


J042450.68+455330.2 


M3 


M4 


- 


M3.5+0.5 


21.6 


122+30 


J183034.15+003800.6 


M6 


M8 


M5.5 


M6.5+1.0 


185.8 


161+50 


J183753.25+001849.2 


M6 


M7 


M5 


M6+1.0 


16.3 


94+30 


J192656.04+21 1438.0 


MO 


MO 


M<2.5 


M0±0.5 


35.6 


52+5 


J202050.40+394243.6 


M2 


M2 


M3 


M2±0.5 


83.5 


140+15 


J202434.30+422126.4 


M5 


M6 


M5.5 


M5.5+0.5 


23.5 


32+3 


J202437.28+385806.9 


M5 


M6 


M4.5 


M5+0.5 


93.2 


138+30 


J202455.53+424504.0 


M5 


M5 


M5 


M5+0.5 


88.6 


127+30 


J202759.82+390418.2 


M3 


M4 


- 


M3.5+0.5 


106.2 


143+10 


J204218.54+395723.8 


M2 


M2 


- 


M2+0.5 


14.5 


19+2 


J204350.68+400108.8 


Ml 


M2 


- 


Ml.5+0.5 


44.2 


71+8 


J204704.82+434911.4 


Ml 


M2 


M2.5 


M2+0.5 


181.0 


245+50 


J205613.08+443424.2 


M6 


M6 


M5 


M6±0.5 


210.2 


261+80 


J205701.63+434138.7 


M7 


M8 


M5.5 


M7+1.0 


40.8 


48.5+20 


J205702.69+434143.7 


M5 


M6 


M5.5 


M5.5+0.5 


117.3 


251+50 


J2 10404.87+535 124.4 


M6 


M6 


M6 


M6+0.5 


282.4 


245+90 




M6 


M6 


M5.5 


M6+0.5 


288.6 


580+100 




M6 


M7 


M5.5 


M6+0.5 


456.1 


570+200 


J213528.41+575823.0 


M4 


M4 


M4.5 


M4+0.5 


126.2 


202+40 


J213545. 87+573640.1 


M3 


M4 


M4 


M4+0.5 


189.8 


221+80 


J213938. 83+575451.4 


M4 


M4 


M5 


M4.5±0.5 


105.4 


198+30 




M4 


M4 


M4.5 


M4±0.5 


124.0 


166+40 


J214547.73+564845.7 


M4 


M6 


M5 


M5+1.0 


147.8 


131+55 


J2 14625.99+572829.0 


M2 


M3 


M2.5 


M2.5+0.5 


243.8 


360+60 




M4 


M5 


M4.5 


M4.5+0.5 


388.3 


541+180 


J222025.22+605423.7 


M4 


M4 


M4.5 


M4+0.5 


233.8 


370+100 


J224830.78+61 1417.9 


M5 


M5 


M4.5 


M5+0.5 


137.7 


204+70 




M5 


M5 


M5 


M5+0.5 


143.7 


190+30 


J22533 1.29+623543.3 


M6 


M7 


M5 


M6+1.0 


194.6 


350+100 


J225923.71+614138.6 


M5 


M5 


M5 


M5+0.5 


250.1 


540+50 



In Fig. [T] we display an example of the application of the;^^ 
method for the spectral type determination in one of our targets. 
The object spectrum (black line) is compared with the template 
that yields the minimum value (red line), and with the tem- 
plates that have adjacent earlier or later spectral subclass (blue 
and green lines, respectively). In general the method works fine, 
although there were some objects that did not fit properly to any 
spectral template. This method did not provide trustworthy re- 
sults for noisy or highly reddened spectra because the calcula- 
tion of is based on the difference point to point between two 
spectra. 



4.2. Equivalent widths 
4.2.1. Spectral indices 

The spectral indices PC3, PC6, TiOl and Ti02 of the objects 
with emission in Ha have been determined using the IRAF task 
stands and the results (see Table [1] only available i n electronic 
form) have been compared with those of iMartm et al. ( 1999) to 
determine the spectral type of our candidates. For 7 objects we 
could not find a reasonable agreement between the different in- 



dices (within 2 spectral subclasses), and hence this method was 
not used for them. 

Hammer provides measurements of the equivalent width 
of Ha simultaneously with the spectral type determination. 
However, it does not provide reliable errors for our spectra be- 
cause Hammer calculates them by propagating the error esti- 
mates provided by the SDSS spectral reduction pipeline. 

We derived the Ha emission equivalent width of our targets 
by direct integration of the line profile using the IRAF task splot. 
The uncertainty was estimated via measurements performed for 
extreme choices of continuum level and line integration wave- 
length range made by visual inspection of the spectra. 

We have tried to measure the Nal doublet equivalent width at 
818.3 nm and 819.5 nm on all of our objects as low surface grav- 
ity indicator, but this part of the spectrum is affected by residual 
fringing and we can not give a reliable estimation of the equiv- 
alent width. To investigate the surface gravity and for compar- 
ison purposes, simulations of how the residual fringing affects 
the doublet on the reference objects have been done. We have 
degraded the resolution of the reference objects to the resolution 
of our targets and introduced different residual fringing accord- 
ing to the residuals in our objects. For statistical and comparison 
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Fig. 2. Comparison of the spectra of IPHAS 
J214625.99+572829.0 observed on two different nights. 
There is significant variabiHty in the continuum and the Ha 
emission. The spectra are normaUsed at 8123A. 



purposes with our targets also measurements of Nal doublet in- 
dex with defined blue, central and red bandpasses (8 130-8 171 A, 
8172-8207 A, 8235-8265A respectively due to low resolution of 
the spectra and visual inspection) have been used for quantita- 
tive estimation of the variation with the residual fringing on the 
objects. The results are shown and discussed in section [53] 

5. Results 

5.1. Definition of the sample of M-type Ha emitters 

9 of our 42 targets with confirmed Ho- emission have very red 
continuum slopes on the red part of the spectrum, indicative of 
high reddening and are classified as early-M or late-K objects. 
They do not have late-M spectral type and thus are not consid- 
ered in the rest of this paper. Our results are focused on the 33 
M-type Ha emission objects that have negligible reddening ac- 
cording to our visual inspection of the comparison with refer- 
ence dwarfs of known spectral types. 

In Table |2] we give our results obtained from the 3 different 
methods of determination from the spectral type (SpT) and the 
measurements of equivalent widths (W) of Ha (using Hammer 
and IRAF) and the Na I doublet. 

6 of our targets were observed more than once to check for 
variability. They have more than one data line in Table |2] 
We find that observations of the same objects obtained at 
different epochs and with different instruments give consistent 
spectral types within the error bars except for object IPHAS 
J214625.99+572829.0 which shows significant variability 
(Fig |2l4b .This object also displays strong and variable Ho- 
emission. Two more objects (IPHAS J210404.87+535124.4, 
J001655. 91+654732.8) show Ha variabiHty larger than the 
error bars but no significant change in spectral type. The three 
remaining targets observed more than once do not have Ha 
variability larger than the error bars. These examples indicate 
that the observed Ha emission is not due to occasional flares, 
because the strong Ha emission is observed at several epochs. 
Emission line variability is an additional indicator for youth and 
has been found previously for young, accreting brown dwarfs 
dScholz & Javawardhanal[2O06h . The variability might be due to 
co-rotating hot spots, changes in the accretion flow geometry 
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Fig. 3. Comparison of spectral type determinations using differ- 
ent methods. 



or changes in the mass accretion rate onto the central object 
(lBouvieretal.|[T995,) . 

In Fig.[3]the results of the three different methods of spectral 
type determination are compared. Hammer tends to give earlier 
types than the method. The discrepancies between these two 
methods is in general within one spectral subclass, but there are 
a few examples of larger differences. The spectral types obtained 
from the spectral indices are sometimes intermediate between 
Hammer and the method, but for the latest subclasses in our 
sample they tend to agree more with Hammer than with the;^'^. 
We have adopted spectral subclasses by computing the average 
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SpT 



Fig. 4. Ha equivalent width against adopted spectral type for our 
objects. The dashed line denotes the dividing line between chro- 
mospheric activity and disk accretion. Our objects are clearly 
above the dashed line and hence they are likely undergoing mass 
accretion. 
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Fig. 5. Effects of the spectral resolution and residual fringing 
on our M6 template (SDSSJ150026.34-003928.0). Dotted line: 
original spectrum of the template; solid line: the same spectrum 
degraded to our spectral resolution; dashed line: the same spec- 
trum degraded to our spectral resolution, and with residual fring- 
ing. For comparison we also show the spectrum of a typical M6 
object from our sample (IPHAS J210404.87+535124.4). 



of the 3 methods, and we have rounded to the nearest subclass 
in steps of 0.5 subclasses. 

With respect to the measurements of equivalent width 
of Ha it is observed that the values derived with Hammer 
are systematically lower than the ones obtained by manual 
integration of the line profile. We consider that the latter ones 
are more reliable than the Hammer results because they rely 
on careful visual inspection of the continuum and line profile 
on each spectrum. Hence, we adopt the manual values of Ha 
equivalent widths and error bars for the subsequent discussion. 



5.2. D\sk accretion 

In order to discern whether the observed Ha emission is likely 
due to chromospheric activity or to disk accretion, we plot our 
Ha equivalent widths versus the adopted spectral type, and 
we compare with the empirical upper limit boundary of chro- 
mosph eric activity derived in Barrado v Navascues & Martini 
(I2003h . The results are shown in Fig. |4l The 6 objects observed 
more than once are marked with red squares and the different 
measurements are connected with red lines. As mentioned in the 
previous section, the persistence of the strong Ha emission indi- 
cates that it is not due to sporadic flares. All the objects occupy 
the region of the diagram above the chromospheric activity satu- 
ration line, and hence they are likely accreting matter from disks. 

5.3. Surface gravity 

The Nal subordinate doublet at 818.3nm and 8 19.5nm is a good 
gravity indicator for late-M spectral types (iKirkpatrick et al.l 
ll99U]Martinet al.lll996t iBriceno et al.1ll998h . This doublet be- 
comes weaker as the surface gravity decreases. We have checked 
whether the doublet can be identified on the reference objects 
despite of the fringing and low resolution of our targets. In Fig. 
IS we compare the SDSSJ150026.34-003928.0 reference dwarf 
(solid and dotted top lines and dashed bottom line) with one of 




Fig. 6. Dependence of the W(Na I) on the residual fringing and 
the spectral type. Filled circles denote the IPHAS objects and 
crosses denote measurements with different fringing residuals 
for the M objects used as reference to classify the IPHAS ob- 
jects. 



our M6 classified objects IPHAS J210404.87-H535124.4 (solid 
bottom line) and the doublet can be clearly be identified despite 
the residual fringing and low resolution. In Fig.|6]we show the 
results of the measurements described in |4.2| We find that for M3 
type and later the doublet can be clearly identified but not for ear- 
lier types since the doublet is weaker for this spectral types and 
is embedded on the residual fringing. According these results we 
can conclude that in all of our objects classified as M3-4 the Nal 
doublet is weaker in our Ha emission objects than in field dwarfs 
of similar spectral class and much weaker for objects later than 
M4. This is consistent with our previous conclusion that most of 
our objects are very young and actively accreting. The half life- 
time for disk accret ion in VLM ob j ects h as been estimated to be 
less than 5 Myr bv lMohantv et all ( l2005h . According to this, the 
results suggest that our sample is <5 Myr. 
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Fig. 7. Colour-magnitude diagrams of the candidates. Open circles represent the objects observed spectroscopically; plus signs 
identify the observed objects that show Ha emission and are classified as M4 or later. 
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Fig. 8. Colour-colour diagrams. Open circles represent all objects observed spectroscopically; open triangles represent the objects 
that show Via emission and are classified as earlier than M4, and plus signs represent the objects that show Ha emission and are 
classified as M4 or later. 
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Fig. 9. Comparison of the results for the IPHAS objects (filled circles) and the characterisation for M dwarfs of lHawlev et al.l (|2002|) 
(open squai^es). 
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Table 3. Photometric data of the IPHAS objects. 



IPHAS name/position 


Observation 


IPHAS magnitudes 






2MASS magnitudes 




J[RA(2000)+Dec.(2000)] 


date 


r' 


r'-i' 


Hff 


J 


H 


K 


J00I649.57+654241.8 


2004-08-07 


20.308+0.045 


2.015+0.054 


2.107+0.050 


15.005+0.053 


13.758+0.052 


12.981+0.041 


J001655.91+654732.8 


2004-08-07 


18.803+0.014 


2.209+0.017 


1.593+0.018 


13.645+0.026 


12.506+0.023 


11.993+0.022 


JOl 1443.00+620820.9 


2003-10-10 


20.033+0.148 


1.756+0.155 


2.115+0.153 


15.487+0.064 


14.198+0.070 


13.411+0.042 


J012348.67+614931.8 


2003-10-12 


19.716+0.081 


1.722+0.086 


2.102+0.083 


15.636+0.061 


14.541+0.060 


14.039+0.053 


J013720.01+645957.7 


2003-10-14 


18.001+0.028 


1.721+0.030 


1.318+0.031 


13.621+0.028 


12.687+0.036 


12.083±0.023 


J023616.00+615609.6 


2003-11-10 


19.727+0.087 


1.280+0.096 


2.041+0.090 


15.984+0.099 


15.022+0.103 


I4.285±0.076 


J035449.17+530903.3 


2004-12-01 


20.249+0.059 


1.961+0.071 


1.906+0.072 


15.719+0.087 


14.719+0.088 


I4.055±0.072 


J035823.95+522312.6 


2004-12-29 


19.016+0.038 


1.283+0.057 


1.965+0.042 


15.745+0.061 


14.826+0.061 


I4.298±0.064 


J042450.68+455330.2 


2004-10-22 


18.619+0.017 


1.732+0.021 


1. 662+0.021 


14.116+0.030 


12.76 +0.026 


I1.849±0.022 


J183034.15+003800.6 


2004-06-08 


18.368+0.010 


1.907+0.013 


1.983+0.012 


13.633+0.028 


12.524+0.031 


11.937±0.024 


J183753.25+00I849.2 


2004-06-09 


19.134+0.018 


2.62I+0.02I 


1.566+0.024 


13.05 ±0.026 


11.837+0.021 


11.16+0.022 


J192656.04+21 1438.0 


2004-08-06 


16.667+0.005 


1. 195+0.008 


1.313+0.007 


12.347+0.022 


11.39 +0.021 


10.742+0.019 


J202050.40+394243.6 


2007-06-21 


19.954+0.037 


I.I 10+0.056 


1.757+0.049 


15.83 ±0.076 


14.569+0.051 


13.611+0.050 


J202434.30+422 126.4 


2003-08-10 


18.304+0.040 


1.912+0.042 


1. 170+0.046 


13.335+0.024 


12.249+0.027 


11.501+0.023 


J202437.28+385806.9 


2006-11-04 


19.879+0.070 


1. 921+0.078 


1.664+0.080 


15.337+0.055 


14.293+0.052 


13.685+0.055 


J202455.53+424504.0 


2004-08-08 


18.902+0.014 


1.910+0.018 


1.525+0.020 


14.425+0.037 


13.548+0.038 


12.883+0.030 


J202759.82+390418.2 


2006-11-24 


17.947+0.013 


1.450+0.017 


1.820+0.015 


13.729+0.030 


12.603+0.026 


11.897+0.022 


J204218.54+395723.8 


2007-06-29 


17.712+0.009 


1.468+0.012 


0.729+0.015 


14.152+0.028 


13.099+0.029 


12.587+0.026 


J204350.68+400108.8 


2007-06-29 


17.273+0.007 


1.270+0.010 


1.223+0.010 


14.035+0.031 


13.024+0.023 


12.44 +0.026 


J204704.82+434911.4 


2003-10-12 


18.758+0.032 


1.232+0.037 


2.049+0.034 


15.336+0.054 


14.343+0.051 


13.764+0.043 


J205613.08+443424.2 


2003-11-10 


20.250+0.113 


1.958+0.121 


2.267+0.120 


14.678+0.043 


13.895+0.056 


13.387+0.045 


J205701.63+434138.7 


2007-06-23 


18.764+0.013 


2.269+0.016 


1.373+0.019 


13.952+0.029 


13.142+0.032 


12.649+0.029 


J205702.69+434143.6 


2007-06-23 


19.871+0.029 


2.252+0.034 


1.868+0.036 


14.348+0.055 


13.382+0.041 


12.695+0.039 


J210404.87+535124.4 


2006-10-02 


18.299+0.015 


1.464+0.020 


2.507+0.016 


13.432+0.025 


12.401+0.025 


11.613+0.025 


J213528.41+575823.0 


2004-08-05 


18.640+0.015 


1.360+0.025 


1.980+0.019 


15.15 ±0.059 


14.212+0.056 


13.598+0.044 


J213545. 87+573640.1 


2004-08-05 


18.526+0.014 


1.509+0.021 


1.996+0.017 


14.551+0.037 


13.676+0.034 


13.143+0.033 


J213938. 83+575451.4 


2004-08-27 


19.637+0.044 


1.658+0.053 


2.083+0.049 


14.757+0.039 


13.621+0.040 


13.032+0.033 


J214547.73+564845.7 


2004-08-28 


18.942+0.035 


2.058+0.039 


1.740+0.040 


14.333+0.032 


13.399+0.036 


12.822+0.025 


J214625.99+572829.0 


2004-08-29 


18.592+0.036 


1.377+0.046 


2.315+0.038 


14.069+0.034 


13.064+0.033 


12.301+0.026 


J222025.22+605423.7 


2003-11-03 


19.167+0.051 


1.459+0.057 


2.130+0.053 


15.32 ±0.053 


14.413+0.051 


13.911+0.052 


J224830.78+611417.9 


2004-08-05 


19.650+0.034 


2.069+0.040 


1.988+0.040 


14.959+0.043 


13.967+0.040 


13.458+0.036 


J22533 1.29+623543. 3 


2004-08-21 


20.554+0.051 


2.153+0.062 


1.997+0.060 


15.393+0.059 


14.17+0.046 


13.471+0.035 


J225923.71+614138.6 


2007-06-27 


19.228+0.016 


1.493+0.023 


2.219+0.019 


15.34I+0.07I 


14.192+0.057 


13.487+0.051 



5.4. Photometric properties 

The main emphasis of our analysis has been on the spectroscopic 
properties of our sample. However, it is also interesting to study 
the photometric properties. The photometry according to best 
new IPHAS photometry of our 33 main objects is given in Table 
[3] colour-magnitude (Fig.|2l) and colour-colour (Fig.[8]l diagrams 
are used to check if our sample of 33 interesting objects stand out 
from the crowd. 

In Fig|7]we have represented with a plus sign the objects classi- 
fied as M4 or later It is observed that they tend to populate the 
fainter part of the diagram. All of them have magnitudes of r' > 
16.5, J >12, i' > 15.5andH > 11. 

In Fig. [8] we plot also the objects that show emission in Ho' but 
are classified as earlier than M4 as open triangles and those clas- 
sified as M4 or later with plus sign. All the accretors show red 
colours with H-K > 0.5, and most of them i'-J < 4. We find 
that the r' -i' and r' -H colours tend to the blue, which could be 
related to the presence of optical veiling in the r' band. On the 
other hand, the red H-K colours could be ascribed to warm dust 
emission. 

We have compared the colour-spectral type relation for our 
targets classified as M4 or later with those in the study of 
iHawlev et alj (|2002|) in which the spectroscopic characterisation 
of SDSS M dwarfs was done (see Fig.|9l). In general our objects 
follow the relation observed in the SDSS objects, but there are 
some significant discrepancies. In the r'-i' vs. SpT diagram, 
it is observed a tendency to the blue of some objects compared 



with SDSS locus. On the other hand, for the colour i'-J, it is ob- 
served a tendency to the red for our sample. These discrepancies 
could be due to the effects of accretion and low surface gravity 
on the spectral energy distribution of VLM objects. 
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Table 1. Spectral indices obtained for the IPHAS candidates. 



Object 


PC3 


PC6 


TiOl 


Ti02 


J001649.57+654241.8 


1.20 


9.50 


1.52 


1.22 


J001655.91+654732.8 


1.23 


8.08 


1.55 


1.19 




1.16 


7.68 


1.44 


1.10 


J011443.00+620820.9 


1.19 


6.47 


1.15 


1.02 


J012348.67+614931.8 


1.13 


4.01 


1.21 


1.05 


J013720.01+645957.7 


1.11 


4.14 


1.30 


1.12 


J023616.00+615609.6 


0.99 


2.54 


1.17 


0.94 


J035449. 17+530903.3 


1.06 


3.59 


1.16 


1.03 


J035823.95+522312.6 


1.00 


2.01 


1.03 


0.88 


J042450.68+455330.2 


1.21 


3.48 


1.14 


0.93 


J183034. 15+003800.6 


1.35 


7.46 


1.22 


0.92 


J183753.25+001849.2 


1.35 


5.77 


1.32 


1.17 


J192656.04+21 1438.0 


1.05 


1.76 


1.03 


0.96 


J202050.40+394243.6 


1.11 


2.42 


0.99 


0.98 


J202434.30+422 126.4 


1.06 


6.91 


1.66 


1.29 


J202437.28+385806.9 


1.12 


4.74 


1.37 


1.18 


J202455.53+424504.0 


1.12 


4.94 


1.56 


1.31 


J202759.82+390418.2 


1.12 


3.00 


1.06 


0.95 


J204218. 54+395723. 8 


1.00 


2.77 


1.04 


0.99 


J204350.68+400108.8 


0.92 


1.78 


1.05 


1.00 


J204704.82+434911.4 


1.08 


2.12 


1.05 


0.97 


J205613.08+443424.2 


1.28 


7.10 


1.28 


1.36 


J205701.63+434138.7 


1.28 


8.90 


1.75 


1.46 


J205702.69+434143.7 


1.09 


7.43 


1.72 


1.33 


J210404.87+535124.4 


1.43 


7.95 


1.26 


1.14 




1.17 


7.09 


1.45 


1.22 




1.30 


7.59 


1.47 


1.27 


J213528.41+575823.0 


1.10 


3.58 


1.15 


1.05 


J213545.87+573640.1 


1.09 


2.97 


1.23 


1.01 


J213938.83+575451.4 


1.09 


4.01 


1.23 


1.03 




1.13 


3.62 


1.23 


1.06 


J214547.73+564845.7 


1.29 


6.43 


1.42 


1.22 


J214625.99+572829.0 


1.05 


3.24 


1.14 


1.08 




1.18 


5.15 


1.27 


1.20 


J222025.22+605423.7 


1.11 


3.11 


1.24 


1.03 


J224830.78+61 1417.9 


1.15 


5.36 


1.50 


1.18 




1.14 


5.70 


1.45 


1.13 


J22533 1.29+623543.3 


1.26 


8.62 


1.36 


1.06 


J225923.71+614138.6 


1.23 


4.75 


1.05 


0.91 



